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Abstract 


Real  radar  data  has  been  analysed  using  the  Fourier  transform  method  and  the  type-IH 
Wigner  distribution.  The  results  show  that  whenever  the  target  was  detectable  by  the 
Fourier  transform  method,  the  target  was  also  detectable  by  the  smoothed  type-HI 
Wigner  distribution  method.  In  the  other  trials,  the  target  was  not  detectable  by  the 
Fourier  transform  method  but  it  was  detectable  by  the  smoothed  type-HI  Wigner 
distribution  method  to  varying  degrees  of  success.  Based  on  the  analysis  of  real  radar 
data,  we  conclude  that  the  smoothed  type-HI  Wigner  distribution  provides  a  promising 
method  of  detecting  accelerating  targets.  However,  more  work  needs  to  be  done  to  find 
an  optimum  smoothing  method.  It  may  turn  out  that  different  smoothing  methods  have 
to  be  used  depending  on  the  acceleration  and  the  closeness  of  the  target  to  the  clutter 
region.  Another  important  contribution  of  the  present  work  is  the  use  of  the  type-HI 
Wigner  distribution  rather  than  the  type-I  Wigner  distribution  which  has  been  used  by 
many  other  researchers.  When  the  type-I  Wigner  distribution  is  used,  the  range  of 
unambiguously  measurable  normalized  velocities  is  tt.  Moreover,  targets  that  are  tt 
radians  away  from  the  clutter  region  in  the  spectral  domain  will  be  masked  by  the 
clutter  and  cannot  be  detected.  The  type-IH  Wigner  distribution  helps  us  to  overcome 
these  problems.  When  the  type-IH  Wigner  distribution  is  used,  the  range  of 
unambiguously  measurable  normalized  velocities  is  27r,  and  target  that  are  tt  radians 
away  from  the  clutter  region  can  be  detected. 
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Resume 


Des  donnas  radar  reelles  ont  dte  analysdes  au  moyen  de  la  mdthode  de  la  transformde 
de  Fourier  et  de  la  distribution  de  Wigner  de  type-III.  Les  rdsultats  montrent  que, 
lorsque  la  cible  etait  detectable  par  la  mdthode  de  la  transformde  de  Fourier,  elle  I’etait 
aussi  par  la  mdthode  de  distribution  de  Wigner  de  type-III  lissde.  Dans  les  autres  essais, 
la  cible  n  eetait  pas  detectable  par  la  mdthode  de  la  transformee  de  Fourier,  mais  elle 
I’dtait  par  la  methode  de  distribution  de  Wigner  de  type-111  lissfe  avec  divers  degres  de 
succds.  L’ analyse  des  donndes  radar  reelles  permet  de  conclure  que  la  distribution  de 
Wigner  de  type-III  lissde  constitue  une  mdthode  prometteuse  de  detection  des  cibles  en 
acceleration.  II  faudra  cependant  mener  d’autres  travaux  pour  trouver  une  methode  de 
lissage  optimale.  II  se  peut  qu’il  faille  utiliser  differentes  methodes  de  lissage  selon 
I’accdleration  de  la  cible  et  sa  proximitd  de  la  rdgion  de  clutter.  Une  autre  contribution 
importante  des  travaux  en  cours  est  I’utilisation  de  la  distribution  de  Wigner  de  type-III 
de  preference  a  la  distribution  de  Wigner  de  type-1,  dont  de  nombreux  autres  chercheurs 
se  sont  servi.  Lorsque  la  distribution  de  Wigner  de  type-I  est  utilisee,  la  fourchette  des 
vitesses  normalisdes  mesurables  non  ambigues  est  tt.  En  outre,  les  cibles  h  tt  radians  de 
la  region  de  clutter  dans  le  domaine  spectral  seront  masquees  par  le  clutter,  d’ou 
impossibilite  de  les  detecter.  La  distribution  de  Wigner  de  type-III  permet  de  surmonter 
ces  problemes.  Lorsque  la  distribution  de  Wigner  de  type-III  est  utilisee,  la  fourchette 
des  vitesses  normalises  mesurables  non  ambigugs  est  2n,  et  les  cibles  h  tt  radians  de  la 
region  de  clutter  peuvent  etre  detectfes. 
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Executive  summary 


One  of  the  central  problems  in  High  Frequency  (HF)  radar  data  is  the  analysis  of  a  time 
series.  The  Fourier  transform  method,  or  Doppler  processing  method,  has  been 
generally  used  in  HF  radar  to  detect  targets  that  are  moving  with  constant  radial 
acceleration.  Examples  of  accelerating  targets  are  manoeuvring  aircrafts  and  missiles. 
Thayaparan  and  Yasotharan  [1]  show  that  there  are  limitations  and  shortcomings  in  the 
Fourier  transform  method  to  detect  accelerating  targets  because  of  the  phenomenon 
known  as  Doppler  smearing.  They  show  that  when  the  target  is  constantly  accelerating, 
the  Fourier  method  may  still  be  used  to  detect  the  target  and  estimate  its  median 
velocity,  provided  the  acceleration  is  small  enough  in  the  sense  to  be  described  in  this 
report.  It  is  shown  that  for  a  given  acceleration,  the  number  of  pulses  cannot  be 
increased  indefinitely  without  resulting  in  catastrophic  failure  of  the  method. 
Conversely,  for  a  given  number  of  pulses,  the  acceleration  cannot  be  arbitrarily  large 
without  resulting  in  catastrophic  failure  of  the  method.  Thus  the  number  of  pulses  and 
the  acceleration  have  to  be  matched  to  achieve  optimum  performance  [1]. 

Consequently,  for  the  interpretation  of  radar  data  in  terms  of  a  changing  frequency 
content,  we  need  a  representation  of  our  data  as  a  function  of  both  time  and  frequency. 
The  purpose  of  this  report  is  therefore  to  stress  the  importance  of  alternative  methods 
which  have  received  little  attention  in  the  past,  namely,  the  joint  time-frequency 
representation  of  signals.  In  this  report,  we  choose  the  Wigner  distribution,  together 
with  a  ambiguity  function,  as  the  basic  tools,  since  both  possess  some  important 
properties  which  make  them  very  attractive  for  time-frequency  signal  analysis  [2]. 

Real  HF  radar  data  has  been  analysed  using  the  Fourier  transform  method  and  the 
type-in  Wigner  distribution.  The  results  show  that  whenever  the  target  was  detectable 
by  the  Fourier  transform  method,  the  target  was  detectable  also  by  the  smoothed 
type-in  Wigner  distribution  method.  In  the  other  trials  the  target  was  not  detectable  by 
the  Fourier  transform  method  but  the  target  was  detectable  by  the  smoothed  type-Hl 
Wigner  distribution  method  to  varying  degrees  of  success.  Based  on  the  analysis  of  real 
radar  data,  we  conclude  that  the  smoothed  type-HI  Wigner  distribution  provides  a 
promising  method  of  detecting  accelerating  targets.  However,  more  work  has  to  be 
done  to  find  an  optimum  smoothing  method.  It  may  turn  out  that  different  smoothing 
methods  have  to  be  used  depending  on  the  acceleration  and  the  closeness  of  the  target 
to  the  clutter  region.  Another  important  contribution  of  the  present  work  is  the  use  of 
the  type-in  Wigner  distribution  rather  than  the  type-I  Wigner  distribution  which  has 
been  used  by  many  other  researchers.  When  the  type-I  Wigner  distribution  is  used,  the 
range  of  unambiguously  measurable  normahzed  velocities  is  tt.  Moreover,  targets  that 
are  tt  radians  away  from  the  clutter  region  in  the  spectral  domain  be  masked  by  the 
clutter  and  cannot  be  detected.  The  type-IH  Wigner  distribution  helps  us  to  overcome 
these  problems.  When  the  type-IH  Wigner  distribution  is  used,  the  range  of 
unambiguously  measurable  normalized  velocities  is  27r,  and  target  that  are  tt  radians 
away  from  the  clutter  region  can  be  detected. 
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Sommaire 


Un  des  principaux  problemes  eprouve  avec  les  donnees  de  radar  haute  frequence  (HF) 
est  I’analyse  d’une  sdrie  chronologique.  La  mdthode  d  transformde  de  Fourier,  ou  la 
methode  de  traitement  Doppler,  a  generalement  ete  utilisee  dans  le  radar  HF  pour  la 
detection  des  cibles  d  acceleration  radiale  constante.  Les  cibles  presentant  une 
acceleration  sont,  par  exemple,  les  aeronefs  qui  effectuent  des  manuvres  et  les  missiles. 
Thayaparan  et  Yasotharan  [1]  montrent  que  la  mdthode  d  transformde  de  Fourier 
comporte  des  limites  et  des  lacunes  pour  la  detection  des  cibles  a  acceleration,  en  raison 
du  phdnomene  appele  la  persistance  Doppler.  Ils  montrent  que  dans  le  cas  d’une  cible  d 
accdldration  constante,  la  methode  de  Fourier  peut  encore  6tre  appliquee  pour  detecter 
la  cible  et  estimer  sa  vitesse  mediane,  d  condition  que  1’ acceleration  soil  assez  faible 
dans  le  contexte  du  present  rapport.  On  montre  que  pour  une  acceleration  donnee,  le 
nombre  d’impulsions  ne  peut  pas  etre  augmentd  indefiniment  sans  qu’il  en  rdsulte  un 
echec  catastrophique  de  la  methode.  Inversement,  pour  un  nombre  donne  d’impulsions, 
I’acceieration  ne  peut  pas  augmenter  arbitrairement  sans  qu’il  en  rdsulte  aussi  un  dchec 
catastrophique  de  la  methode.  Ainsi,  le  nombre  d’impulsions  et  I’acceieration  doivent 
etre  adaptds  pour  peimettre  des  performances  optimales  [1]. 

Par  consequent,  pour  interpreter  les  donnees  radar  du  point  de  vue  d’un  contenu  en 
frequences  variable,  nous  devons  les  reprdsenter  aussi  bien  en  fonction  du  temps  qu’en 
fonction  de  la  frequence.  L’objet  du  present  rapport  est  done  de  souligner  I’importance 
des  rndthodes  de  rechange  auxquelles  on  a  accorde  peu  d’ attention  par  le  passe, 
notamment  la  representation  tenq)s-frequence  combinee  des  signaux.  Dans  le  rapport, 
nous  choisissons  la  distribution  de  Wigner  avec  une  fonction  d’ambiguite  comme  outil 
de  base,  etant  donne  que  ces  deux  elements  offrent  certaines  proprietes  importantes  qui 
presentent  un  grand  intdret  pour  1’ analyse  temps-frequence  des  signauxa[2]. 

Des  donnees  radar  HF  rdelles  ont  dtd  analysdes  au  moyen  de  la  methode  de  la 
transformee  de  Fourier  et  de  la  distribution  de  Wigner  de  type-in.  Les  rdsultats 
montrent  que,  lorsque  la  cible  dtait  detectable  par  la  rndthode  de  la  transformee  de 
Fourier,  elle  I’etait  aussi  par  la  methode  de  distribution  de  \Ntigner  de  type-IH  Ussde. 
Lors  des  autres  essais,  la  cible  n’dtait  pas  detectable  par  la  methode  de  la  transformee 
de  Fourier,  mais  elle  I’etait  par  la  rndthode  de  distribution  de  Wigner  de  type-III  Ussde 
avec  divers  degrds  de  succ^s.  L’ analyse  des  donndes  radar  rdelles  permet  de  conclure 
que  la  distribution  de  Wigner  de  type-HI  lissde  constitue  une  rndthode  prometteuse  de 
detection  des  cibles  en  acceleration.  II  faudra  cependant  mener  d’ autres  travaux  pour 
trouver  une  mtithode  de  lissage  optimale.  II  se  peut  qu’il  faille  utiliser  difftirentes 
methodes  de  hssage  selon  I’acceieration  de  la  cible  et  sa  proximite  de  la  region  de 
clutter.  Une  autre  contribution  importante  des  travaux  en  cours  est  FutiUsation  de  la 
distribution  de  Wigner  de  type-IH  de  preference  a  la  distribution  de  Wigner  de  type-I, 
dont  de  nombreux  autres  chercheurs  se  sont  servi.  Lorsque  la  distribution  de  Wigner  de 
type-I  est  utihsee,  la  fourchette  des  vitesses  normalisees  mesurables  non  ambigues  est 
TT.  En  outre,  les  cibles  ^  tt  radians  de  la  region  de  clutter  dans  le  domaine  spectral 
seront  masqudes  par  le  clutter,  d’ou  impossibilite  de  les  ddtecter.  La  distribution  de 
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Wigner  de  type-III  permet  de  surmonter  ces  problemes.  Lorsque  la  distribution  de 
Wigner  de  type-III  est  utilisde,  la  fourchette  des  vitesses  normalisees  mesurables  non 
ambigues  est  2tv,  et  les  cibles  a  tt  radians  de  la  rdgion  de  clutter  peuvent  etre  detectees. 
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1 .  Introduction 


The  Fourier  Transform  is  at  the  heart  of  a  wide  range  of  techniques  that  are  used  in  HF 
radar  data  analysis  and  processing.  Mapping  the  data  into  the  temporal  frequency 
domain  is  an  effective  way  of  recording  the  data  such  that  their  global  characteristics 
can  be  assessed.  However,  the  change  of  frequency  content  with  time  is  one  of  the  main 
features  we  observe  in  HF  radar  data.  Because  of  this  change  of  frequency  content  with 
time,  radar  signals  belong  to  the  class  of  non-stationary  signals. 

One  of  the  central  problems  in  High  Frequency  radar  data  is  the  analysis  of  a  time 
series.  The  Fourier  transform  method,  or  Doppler  processing  method,  has  been 
generally  used  in  HF  radar  to  detect  targets  that  are  moving  with  constant  radial 
acceleration.  Examples  of  accelerating  targets  are  manoeuvring  aircraft  and  missiles. 
Thayaparan  and  Yasotharan  [1]  show  that  there  are  limitations  and  shortcomings  in  the 
Fourier  transform  method  to  detect  accelerating  targets  because  of  the  phenomenon 
known  as  Doppler  smearing.  They  show  that  when  the  target  is  constantly  accelerating, 
the  Fourier  method  may  still  be  used  to  detect  the  target  and  estimate  its  median 
velocity,  provided  the  acceleration  is  small  enough  in  the  sense  to  be  described  in  the 
paper.  It  is  shown  that  for  a  given  acceleration,  the  number  of  pulses  cannot  be 
increased  indefinitely  without  resulting  in  catastrophic  failure  of  the  method. 
Conversely,  for  a  given  number  of  pulses,  the  acceleration  cannot  be  arbitrarily  large 
without  resulting  in  catastrophic  failure  of  the  method.  Thus  the  number  of  pulses  and 
the  acceleration  have  to  be  matched  to  achieve  optimum  performance  [1]. 

Consequently,  for  the  interpretation  of  radar  data  in  terms  of  a  changing  frequency 
content,  we  need  a  representation  of  our  data  as  a  function  of  both  time  and  frequency. 
The  purpose  of  this  paper  is  therefore  to  stress  the  importance  of  alternative  methods 
which  have  had  little  attention  in  the  past,  namely,  the  joint  time-frequency 
representation  of  signals.  In  this  report,  we  use  the  new  discrete-time  \Yigner 
time-frequency  distribution.  This  distribution  was  independently  derived  by  the  authors 
Thayananthan  and  Yasotharan  [3].  But  it  had  already  been  derived  by  Chan  [4]  in  an 
effort  to  solve  the  problem  of  aliasing.  Nevertheless,  the  optimality  properties  of  this 
discrete-time  Wigner  distribution  for  signal  detection  were  not  observed  in  [4]. 
Therefore,  in  the  context  of  signal  detection,  this  discrete-time  \Nfigner  distribution 
seems  new. 
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Z  Signal  Model 


In  this  study,  the  radar  is  assumed  to  be  of  the  Pulse  Doppler  type,  that  is,  it  sends  out  a 
uniform  train  of  RF  pulses  and  phase-coherently  receives  their  returns.  It  is  also 
assumes  that  the  receiver  has  pulse  compression  capability  so  that  a  pulse  return  from 
an  isolated  target  can  be  represented  by  a  single  sample  of  the  compressed  pulse. 


2.1  Single-Target  Scenario 

Suppose  the  radar  sends  out  N  pulses,  one  every  T  seconds,  and  there  is  a  target 
moving  with  constant  radial  acceleration.  Assume  that  the  range  walk  is  negligible,  that 
is,  the  change  in  range  during  the  observation  period  of  NT  seconds  is  negligible 
compared  to  the  radar  range  resolution  as  determined  by  the  width  of  the  compressed 
pulse.  Then  the  samples  of  the  N  range-compressed  pulses  taken  at  the  range  of  the 
target  have  the  form 


(1)  r(n)  =  s(n)-|-t;(n),  forn  =  0, 1,2,3,. (AT  -  1), 


where 


(2)  s(n)  = 

is  the  noise-free  signal,  and  v(n)  is  a  sequence  of  independently  and  identically 
distributed  (iid)  samples  of  complex  Gaussian  noise  with  mean  zero  and  variance  cr^. 

The  signal  parameters  a  and  bo  are  the  target  amplitude  and  phase  respectively.  The 
signal  parameters  6i  and  62  are  the  normalized  initial  radial  velocity  and  the  normalized 
radial  acceleration  respectively. 

2.1.1  Normalized  Initial  Velocity  and  Acceleration 

The  normalized  initial  radial  velocity  is  defined  as 


where  u  is  the  initial  radial  velocity  in  meters/sec,  T  is  pulse  repetition 
interval  in  secs,  and  A  is  the  carrier  wavelength  in  meters.  Similarly,  the 
normalized  radial  acceleration  is  defined  as 
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(4) 


where  /  is  the  radial  acceleration  in  meters/sec/sec  towards  the  radar.  It  can  be 
seen  that  both  of  the  above  normalized  quantities  are  non-dimensional. 

To  obtain  these  relations  note  that  the  increment  in  range  d,  measured  from 
the  beginning  of  the  observation  interval,  as  a  function  of  continuous-time  t,  is 
d{t)  =ut  +  and  therefore  the  phase  increment  of  a  pulse  as  a  function 
of  discrete-time,  or  pulse  index,  is  given  by 


(4mT\  l,/47rr2\  3 

In  the  rest  of  the  discussion,  the  term  velocity  and  acceleration  refer  to 
normalized  radial  velocity  and  normalized  radial  acceleration 
respectively. 

2.1 .2  Per-Pulse  Signal-to-Noise  Ratio 

The  Per-Pulse  Signal-to-Noise  Ratio  is  defined  as 

(7)  SNRp„lse=  0)". 

The  performance  measures  of  any  detection  method  will  depend  on  SNRpuise 
and  the  number  of  pulses  N.  The  values  of  the  signal  parameters  bx  and  ^ 
may  also  influence  performance.  We  are  interested  in  detection  methods 
whose  performance  will  increase  as  N  increases  regardless  of  the  values  of  bi 
and  &2- 


2.2  Multiple-Target  Scenario  -  Targets  Separated  in  Range 

When  there  are  multiple  targets,  as  long  as  they  are  separated  from  one  another  in  range 
by  at  least  the  radar  range  resolution  as  determined  by  the  width  of  the  compressed 
pulse,  a  model  of  the  above  form  will  be  valid  for  each  of  them. 

Thus,  if  there  are  two  targets  separated  in  range,  the  samples  of  the  N 
range-compressed  pulses  taken  at  the  corresponding  ranges  of  the  targets  have  the  form 
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(8)  ri(n)  =  si(n)  +  Wi(n),  and 

(9)  r2(n)  =  S2(n)  +  t>2(n) 


forn  =  0,l,2,3,...,(iV-l), 


where  si  (n)  and  S2  (n)  are  the  noise-free  signals,  and  vi  (n)  and  t;2  (n)  are  sequences  of 
independently  and  identically  distributed  samples  of  complex  Gaussian  noise  with 
mean  zero  and  variance 

Because  of  the  independence  of  t;i(n)  and  V2(h),  the  model  decouples  into  two  single 
target  models. 


2.3  Multiple-Target  Scenario  -  Targets  Coincident  in  Range 

When  there  are  multiple  targets  coincident  in  range,  the  samples  of  the  N 
range-compressed  pulses  taken  at  the  range  of  the  targets  have  the  form 


(10)  r(n)  =  s(n)-|-r;(n),  forn  =  0,l,2,3,...,(iV  - 1), 


where  a(n)  is  the  combined  noise-free  signal  due  to  all  targets  and  v(n)  is  a  sequence 
of  independently  and  identically  distributed  samples  of  complex  Gaussian  noise  with 
mean  zero  and  variance  If  there  are  K  targets,  the  combined  noise-free  signal  has 
the  form 


(11) 


K 

s(n)  =  y^a^.e^(^k,o+i>k,in+jbk,2n’‘)^ 
k=l 


where  Ck  and  bkfi  are  the  amplitude  and  phase  respectively  of  the  target,  and  bk^i 
and  bk,2  are  the  normalized  initial  radial  velocity  and  the  normalized  radial  acceleration 
respectively  of  the  target. 
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3.  A  Wigner  Distribution  Formuiation  for  a 

Discrete^Time  Chirp  Signal  _ 

It  is  nearly  20  years  ago  that  the  Wigner  distribution  (also  known  as  Wigner- Ville 
distribution)  was  given  a  prominent  place  in  the  signal  processing  community’s 
spectrum  of  research  interests.  Recognition  was  given  almost  immediately  for  the 
potential  of  the  Wigner  distribution  as  a  tool  for  displaying  and  analyzing  signal 
characteristics  in  the  time-frequency  plane.  This  assessment  was  soundly  supported  by 
the  long  list  of  desirable  mathematical  properties,  such  as  correct  marginal  and  support 
properties,  satisfied  by  the  Wigner  distribution.  Unlike  the  spectrogram,  for  instance, 
the  resolution  of  the  \Nfigner  distribution  is  not  limited  by  the  temporal  or  spectral 
properties  of  a  window,  since  a  window  is  not  required  in  the  definition  of  the  Wigner 
distribution  [3,5]. 

On  the  other  hand,  due  to  its  high  resolution  and  its  quadratic  nature,  the  Wigner 
distribution  locally  exhibits  oscillations  between  positive  and  negative  values  to  an 
extent  that  smoothing,  in  a  certain  sense  in  accordance  with  the  Heisenberg  uncertainty 
principle,  is  necessary  to  make  it  non-negative  everywhere  [6].  As  a  consequence,  a 
pointwise  interpretation  of  the  Wigner  distribution  as  a  true  probability  density  function 
is  not  possible.  This  circumstance  has  kept  the  discussions  in  the  signal  processing 
community  quite  animated,  yielding  all  sorts  of  modified  Wigner  distributions  through 
the  design  of  smoothing  kernels  and  distributions  involving  the  signals  in  a 
non-quadratic  way.  Thus,  to  a  large  extent,  the  interest  and  research  effort  of  the 
community  has  been  directed  at  redefining  and  re-interpreting  the  Wigner  distribution, 
with  the  purpose  of  obtaining  a  clear  and  useful  display  and  analysis  tool. 

Besides  being  a  useful  display  and  analysis  tool,  the  Wigner  distribution  can  be  a  great 
help  in  the  design  of  time-frequency  filtering  methods.  In  this  section  and  the  following 
sections,  we  develop  such  methods  for  the  detection  of  aircraft  from  the  noisy  signals. 
The  Wigner  distribution  was  originally  defined  for  continuous-time  signals.  A 
discrete-time  version  of  it  was  proposed  recently.  Unfortunately,  this  discrete-time 
Wigner  distribution  suffers  from  aliasing  effects,  which  prevent  several  of  the  properties 
of  the  continuous-time  Wigner  distribution  from  carrying  over  straightforwardly.  In  this 
report,  a  discrete-time  Wigner  distribution  which  does  not  suffer  from  aliasing  is 
introduced.  It  is  essentially  an  augmentation  of  the  previous  version,  incorporating  new 
information  about  the  signal  not  contained  in  the  previous  version. 


3.1  Three  Definitions  of  Wigner  Distribution  for  Discrete-Time 
Signais 

Given  a  discrete-time  signal  r(n),  we  consider  three  definitions  of  auto  time-frequency 
distributions  which  have  some  properties  analogous  to  those  of  the  classical  Wigner 
distribution  defined  for  continuous-time  signals  [7]. 
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It  is  straightforward  to  verify  that  these  are  real- valued  functions. 


3.1.1  Type-1  Wigner  Distribution 

The  type-I  Wigner  distribution  W}{n,6)  is  defined  as 


(12)  Wl{n,  r{n  +  jb)r*(n  - 

k 


where  n  is  integer-valued  and  0  is  real-valued.  Note  that  this  is  the  same  as  the 
definition  provided  in  [7]  except  for  the  missing  scaling  factor  2  at  the  front. 

It  is  easy  to  verify  that  the  type-I  Wigner  distribution  is  a  periodic  function  of  9 
with  period  tt.  This  property  are  discussed  more  in  the  following  subsections. 

For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (iV  -  1),  the  type-I  Wigner 
distribution  is  zero  outside  0  <  n  <  (JV  —  1). 

3.1.2  Type-ll  Wigner  Distribution 

The  type-n  Wigner  distribution  6)  is  defined  as 


(13)  W^^{n,d)  =  ;^r(n-hJk-b  l)r*(n-  ifc)e-^(2fc+i)e^ 

k 


where  n  is  integer- valued  and  0  is  real- valued. 

It  is  straight  forward  to  verify  that  the  type-II  Wigner  distribution  is  a  periodic 
function  of  0  with  period  27r.  It  is  also  straight  forward  to  verify  that 


(14)  9  +  n)  =  9). 

These  properties  are  discussed  more  in  the  following  subsections. 

For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (JV  -  1),  the  type-II  Wigner 
distribution  is  zero  outside  0  <  n  <  (JV  —  2). 
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3.1 .3  Type-Ill  Wigner  Distribution 

The  type-m  Wigner  distribution  0)  is  defined  in  terms  of  the  type-I 

and  type-n  Wigner  distributions  as  follows. 


(15) 


(n/2, 0)  for  even  n, 

-l)/2,0)  for  odd  n. 


It  is  straight  forward  to  verify  that  the  type-IH  Wigner  distribution  is  a 
periodic  function  of  0  with  period  27r.  This  property  is  discussed  more  in  the 
following  subsections. 

For  a  signal  r(m)  that  is  zero  outside  0  <  m  <  {N  —  1),  the  the  type- HI 
Wigner  distribution  is  zero  outside  0  <  n  <  2(iV  —  1).  However,  if  we 
consider  even  n  to  correspond  to  integer  values  of  time  and  odd  n  to 
correspond  to  half-integer  values  of  time  then  the  type-DI  Wigner  distribution 
is  zero  outside  the  time  range  0  <  m  <  {N  — 


3.2  The  Range  of  Unambiguously  Measurable  Velocities  of  the 
Type-1  and  Type-lll  Wigner  Distribution  Based  Methods 

As  shown  in  Section  3.1.1,  the  period  of  the  type-I  Wigner  distribution  is  tt.  Therefore 
the  signals 


(16)  aej(&o+!>in+i62n^)^  and 

(j-y)  jjgj(6o+(6i-Hr)n+562n2)^ 

gives  rise  to  the  same  type-I  Wigner  distribution.  In  otherwords,  if  we  detect,  using  the 
type-I  Wigner  distribution,  the  presence  of  a  signal  in  the  time-frequency  plane  with 
intercept  and  slope  parameters  ci  and  ca  respectively,  then  we  cannot  resolve  the 
ambiguity  as  to  whether  the  physical  signal  was 


(18)  ae?(co+«n-|-ic2n2)^  Qj. 

l2g?(co+(ci+ir)n+|c2n^) 

without  external  information. 

In  contrast,  as  shown  in  Section  3.1.3,  the  period  of  the  type>in  Wigner  distribution  is 
27r.  Therefore  the  signals 
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(20) 

(21) 


oej(6o+f>in+i62n=')^  and 
ae^(^ + (&i  +5r)n+ l&2n^ ) 


give  rise  to  different  type-III  Wigner  distributions. 

If  we  detect,  using  the  type-III  Wigner  distribution,  the  presence  of  a  signal  in  the 
time-frequency  plane  with  intercept  and  slope  parameters  ci  and  C2  respectively,  then 
there  is  no  ambiguity  as  to  whether  the  physical  signal  was 


(22)  aej(co+cin+|c2n*)^ 

(23)  ae;(oo+(ci+ir)n+|c2n2) 


However  we  need  external  information  to  resolve  the  ambiguity  as  to  whether  the 
physical  signal  was 


(24) 

(25)  agj(co+(<n+2T)n+ic2n*) 

3.3  The  Wigner  Distributions  of  the  Discrete-Time  Chirp  Signal 

In  this  sub-section  we  derive  the  type-I,  type-II,  and  type-HI  Wigner  distributions  for 
the  discrete-time  chirp  signal 

(26)  s(n)  =  for  n  =  0, 1, . . . ,  (JV  -  1), 

defined  in  Section  2.  These  derivations  have  been  done  in  Appendix  A  assuming  a  =  1. 
The  results  stated  here  are  obtained  by  scaling  those  derived  in  the  appendix  by  o^. 

Note  that  the  Wigner  distributions  do  not  depend  on  5o- 

3.3.1  The  Type-1  Wigner  Distribution  of  s{n) 

The  Type-I  Wigner  Distribution  VT/  (n,  0)  of  s{n)  is  given  by 


Wjin,0)=a^ 

(27) 


2  min(n,  N  -l-n)  +  l 

Bin[(g-(fri+&an))(2min(n.Af-l-n)-H)l 

8inie--(6i+62n)] 


ifd  =  bi+  b2n  mod  n, 
otherwise. 
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for  0  <  n  <  (JV  —  1),  and  it  is  zero  outside  this  range. 

In  a  3-dimensional  plot  [3],  it  has  a  ridge  along  the  lines  given  by 

(28)  6  =  bi+  b^n  mod  tt. 

The  height  of  the  ridge  is  given  by 

(29)  2  min(n,  iV  -  1  -  n)  -|- 1. 

3.3.2  The  Type-ll  Wigner  distribution  of  s(n) 

The  Type-II  Wigner  Distribution  W/^(n,  9)  of  s(n)  is  given  by 


Wl\n,9)=a^^ 

(30) 


2  min(n,  JV  —  2  —  n)  -I-  2 
— (2min(n,  N  —  2  — n)  + 2) 

Bm[(tf— (&l+^fe2+&2w))(2mm(Ti;Ar— 2-n)+2)] 
sin  [0-(bi+^b2-\-h2n^ 


if  5  —  6i  +  mod  27r, 

if  0  ^  n  +  bi  +  hn  mod  27r, 

otherwise. 


for  0  <  n  <  (iV  ~  2),  and  it  is  zero  outside  this  range. 

In  a  3-dimensional  plot  [3],  it  has  a  ridge  along  the  lines  given  by 


(31)  0  =  bi+  ^&2  4-  b2n  mod  27r. 

£t 

The  height  of  the  ridge  is  given  by 


(32)  2  min(n,  JV  -  2  -  n)  -b  2. 

In  a  3-dimensional  plot,  it  also  has  a  valley  along  the  lines  given  by 

(33)  6  =  iT  +  bi+^b2  +  b2n  mod  tt. 

The  depth  of  the  valley  is  given  by 

(34)  2  min(n,  JV  -  2  -  n)  -I-  2. 
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3.3.3 


I 


The  Type-lll  Wigner  distribution  of  s(n) 

The  Type-III  Wigner  Distribution  0)  of  s(n)  is  given  by 


(35) 


2min(f,iV-l-f)  +  l 
(-l)"[2inin(2,JV-l-f)  +  l] 

8m[(tf-(6i+l62n))(2min(§,JV-l-2)+i)l 


if6  =  bi+  ^b2n  mod  27r, 
if  0  =  TT  +  6i  +  mod  27r, 

Otherwise. 


for  0  <  n  <  2(A^  —  1),  and  it  is  zero  outside  this  range. 

In  a  3“dimensional  plot  [3],  it  has  a  ridge  along  the  lines  given  by 


(36)  0  =  bi  +  -b2n  mod  27r. 

The  height  of  the  ridge  is  given  by 


(37)  min(n,  2{N  -  1)  -  n)  +  1. 

In  a  3-dimensional  plot,  it  also  has  an  oscillatory  structure  along  the  lines  give 
by 

(38)  ^  —  TT  +  mod  27r. 

A 

The  period  of  the  oscillation  is  one  unit  and  the  amplitude  of  the  oscillation  is 
given  by 


(39) 


min(n,  2(iV  -  1)  -  n)  +  1. 
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4.  An  Ambiguity  Function  Formulation  for  a 

Discrete-Time  Chirp  Signal _ 

The  inverse  Fourier  transform  of  the  Wigner  distribution  is  called  the  ambiguity 
function  (AF).  The  Fourier  transform  maps  the  Wigner  distribution  auto-components  to 
a  region  centered  on  the  region  of  the  AF  plane,  whereas  it  maps  the  oscillatory  Wigner 
distribution  cross-components  away  from  the  origin  [8]. 

The  fact  that  the  auto-  and  cross-components  are  spatially  separated  in  the  AF  domain 
means  that  if  we  apply  mask  function  to  the  AF,  we  can  suppress  some  of  the 
cross-components.  This  masking  operation  defines  a  new  time-frequency 
representation  (TFR) 


(40)  TFR  =  Fourier  trajasform{AF  •  Kernel} 

with  properties  different  from  the  Wigner  distribution.  The  mask  function  is  called  the 
‘kernel’  of  the  TFR.  Since  there  are  many  possible  2-dimensional  kernel  functions, 
there  exist  many  different  TFRs  for  the  same  signal.  The  class  of  all  TfRs  obtained  in 
this  fashion  is  called  Cohen’s  class.  A  more  detailed  description  of  the  ambiguity 
function  is  given  [8]. 


4.1  Three  Definitions  of  Ambiguity  Function  for  Discrete-Time 
Signais 

Given  a  discrete-time  signal  r(n),  we  consider  three  definitions  of  ambiguity  function 
which  are  related,  via  a  two-dimensional  Fourier  transformation,  to  the  three  definitions 
of  Wigner  distribution  given  in  Section  3.1 . 

4.1.1  Type-I  Ambiguity  Function 

The  type-I  ambiguity  function  A^(A;,  a)  is  defined  as 


(41)  A^(A;,a)  =  ^r-(m  +  k)r*(m  —  A:)e’^’"“, 

m 


where  k  is  integer- valued  and  a  is  real- valued. 

For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (JV  -  1),  the  type-I  ambiguity 
function  is  zero  outside  —floor  (•^^)  <k<  floor 
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4.1.2 


Type-ll  Ambiguity  Function 

The  type-II  ambiguity  function  a)  is  defined  as 


(42)  a)  =  ^r(m  +  k  +  l)r*(m  - 

m 


where  k  is  integer- valued  and  a  is  real- valued. 

For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (AT  -  1),  the  type-II  ambiguity 
function  is  zero  outside  -floor  <  k  <  floor  (^^). 

4.1.3  Type-Ill  Ambiguity  Function 

The  type-lll  ambiguity  function  a)  is  defined  in  terms  of  the  type-I 

and  type-II  ambiguity  functions  as  follows. 


(43) 


A^(k/2,a) 

Al^((k-l)/2,a) 


for  even  k, 
for  odd  k. 


For  a  signal  r(n)  that  is  zero  outside  0  <  n  <  (AT  -  1),  the  type-in  ambiguity 
function  is  zero  outside  —(N  —  1)  <  k  <  (N  —  1). 

4.1 .4  Cross  Ambiguity  Functions  Between  Two  Discrete-Time  Signals 

Given  discrete-time  signals  a;(n)  and  y(n),  we  can  define  cross  ambiguity 
functions  as  straightforward  extensions  of  the  definitions  given  above.  Thus 


(44) 

(45) 

(46) 


Ail^{k,a) 


^a;(m  -I-  k)y*(m  -  k)eP'^°‘, 

m 

^®(m  -F  A:  -f-  l)y*{m  -  , 

m 

/  for  even  k, 

I  1)/2,Q!)  for  odd  k. 


4.2  The  Ambiguity  Functions  of  a  Discrete-Time  Chirp  Signal 

In  this  sub-section  we  derive  the  type-I,  type-II,  and  type-III  ambiguity  functions  of  the 
discrete-time  chirp  signal 
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*(»)  =  {  0 


gj(6in+f62n2)  jfo<n<(JV-l) 

0  otherwise. 


4.2.1  Type-1  Ambiguity  Function 

Recall  that  the  type-I  ambiguity  function  Al{k,  a)  is  defined  as 


Al{k,  a)  =  ^  s{m  +  k)s*{m  — 


where  k  is  integer- valued  and  a  is  real- valued. 

It  is  straight  forward  to  see  that  the  signal  product  term  s(m  +  k)s*(m  —  k)  is 
zero  outside  the  ranges  —floor  <k<  floor  and 

max(A:,  —A;)  <  m  <  iV  —  1  —  max(fe,  —k). 

Thus,  for  —floor  <k<  floor  we  have 


N—l—max{kj—k) 

(49)  =  ^  ^2(b2k+a)m^ 

m=jnsx{kj—k) 

which  is  a  sum  of  a  geometric  series  that  can  be  evaluated  as 


Ai{k,a)  = 


ej2t>ik(jV  _  2max(fc,  —k)) 

ej26ikgi[(62fc+a)(JV-l)] 


smnbok+anis  —zmSiXi 


8iii[(fe2fe+Q!)' 


if  62A;  -h  a  =  0  mod  tt 
otherwise. 


4.2.2  Type-ll  Ambiguity  Function 

Recall  that  the  type-II  ambiguity  function  A\^{k,  a)  is  defined  as 


(51)  Af(A;,a)  =^s(m-hA;-hl)s*(m-ik)e^‘(2"*+i)“ 


where  k  is  integer- valued  and  a  is  real- valued. 
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It  can  be  seen  that  the  signal  product  term  s{m  +  A;  +  l)s*(m  —  k)  is  zero 
outside  the  ranges  —floor  (^)  <k<  floor  and 
max(A:, -{k  +  l))<m<N-2-  max(A;,  -{k  +  1)). 

Thus,  for  -floor  (Jj)  <k<  floor  we  have 


7V-2-inax(A:,-(A:+l)) 

(52)i4f(A;,a)  =  ^  Qj2{b2(k+i}+a)m^ 

m=max(fe,— (fc+l)) 

which  is  a  sum  of  a  geometric  series  that  can  be  evaluated  as 


Ai^iKa) 

(53) 


-  1  -  2max(A:,  -(*  +  1)))  if  b2(A;  + 1)  +  a  =  0  mod  27r, 

ej25i(fc+i)^[(62(fc+i)+a)(Ar-l)]  { sin[(b2(fc+i)+a)(Af-l-2maJc(fc,-(fc+l)))]\  .  . 

*iii[(fc2(fc+i)+a)]  j  otnerwK 


4.2.3  Type-lll  Ambiguity  Function 

Recall  that  the  type-III  ambiguity  function  Al^^{k,a)  is  defined  in  terms  of 
the  type-I  and  type-II  ambiguity  functions  as 


(54) 


Al{k/2,a) 

Ai\ik-l)/2,a) 


for  even 
for  odd 


k, 

k. 


Thus  for  -{N  —  1)  <  k  <  (N  —  1), 


(55) 


e»«>ifc(JV-max(A:,-A;))  if  ^62^  +  0  =  Omod  27r 

ej6ifcej[(|62fc+a)(^-l)]  /  sin[(|62fc+a)(JV-max(fc,-fc))]\  ■ 

\  8m[(i&2fc+Q!)]  J 


4.3  The  Cross  Ambiguity  Functions  of  Two  Discrete-Time  Chirp 
Signals 

In  this  sub-section  we  derive  the  type-I,  type-II,  and  type-III  cross  ambiguity  functions 
of  the  discrete-time  chirp  signals 
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(56) 


(57) 


a:(n)  =  | 
y{n)  =  I 


g}((Hn+io2n^) 

0 


g^(6in+^62n^) 

0 


ifO<n<(Ar-l) 

otherwise. 


ifO<n<  (iV-1) 
otherwise. 


4.3.1  Type-1  Cross  Ambiguity  Function 

Recall  that  the  type-I  cross  ambiguity  function  Al.y{k,  a)  is  defined  as 


(58)  A^yik,  a)  =Y, 

m 


where  k  is  integer- valued  and  a  is  real- valued. 

It  is  straight  forward  to  see  that  the  signal  product  term  x{m  k)y^{m  —  k)  is 
zero  outside  the  ranges  —floor  <  floor  and 

max(fc,  —k)  <  m  <  JV  —  1  —  max(fc,  — fc). 


Defining 


(59) 

di  -\-bi  — 

Cl, 

(60) 

Ci2  -\-b2  — 

C2, 

(61) 

ai  - 

-bi  = 

di, 

(62) 

tt2  ’ 

-62  = 

d2, 

we  can  show 

(63) 


iV — 1 — inax(  fc,— fe) 

E  ^‘(  (<ii +C2  Aj+2Q;)m+ 1  d2m^  ) 

m=:niax(fe,— fe) 


for  —floor  ^  fc  <  floor 


Special  Case  ^2=0 
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(64) 


ejcifc(jy  _  2  inax(A:,  -A:))  if  di  +  C2A:  +  2<i  0  mod  tt, 

^  gJcifcg;[i(<ii+C2A:+2a)(Ar-l)]  ( sinf^ (di+C2k+2a)(N-2 majc(<;,-A;))] "N 
^  y  8in[i(di+c2fc+2Q)]  j  Oi'ierwise. 


Thus  is  concentrated  along  the  line  di  +  C2A;  +  2<i  =  0  and 

symmetrically  spread  about  it  with  respect  to  a.  For  a  fixed  k,  the  amount  of 
spread  with  respect  to  a  can  be  measured  by  1  /{N  -  2  max(fc,  - jfc)). 


The  General  Case 


Although  it  seems  difficult  to  evaluate  A^  y^k,  a)  in  the  general  case,  by 
writing 


Aiy{k,a)  =  eJ(«*+2*''")e^((‘*>+2'^2(V)+':2k+2a)(-*^)) 
iV-l-max(A:,-fc)  . 

^  gf((*+<fc(^)+C2fe+2a))(m-(ii^|a))+ld,(„_(^))a)^ 

Tn=mQx(kj—k) 

(65) 


we  can  see  that  |A^_j,(A:,Q!)|  is  symmetrically  spread  about  the  line 
di  +  d2  +  C2 A;  +  2a  =  0  with  respect  to  a.  For  small  ^2,  | ^x,y{K  o)  | 

is  concentrated  along  the  same  line.  For  a  fixed  k,  the  amount  of  spread  with 
respect  to  a  will  depend  on  both  {N  -2  max(A:,  -A:))  and  d2. 

4.3.2  Type-ll  Cross  Anabiguity  Function 

Recall  that  the  type-II  cross  ambiguity  function  J^^y{k,  a)  is  defined  as 


(66)  £Jy{k,  a)  =  ^  a:(m  +  A;  +  \)y*{m  -  A:)e»(^+*)“, 

m 


where  k  is  integer-valued  and  a  is  real-valued. 

It  is  straight  forward  to  see  that  the  signal  product  term 
x(m  -f  A;  -b  l)y*(Tn  —  A;)  is  zero  outside  the  ranges 
-floor  (y)  <k<  floor  and 

max(A;,  -(A;  -|-l))<m<)V-2  -  max(A:,  -(A:  +  1)). 
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Defining 


(67) 

(68) 

(69) 

(70) 

we  can  show 


0-1+ bi=  Cl, 
02  +  62  =  C2, 
oi  —  61  =  di, 
O2  —  62  ==  ^2, 


AT— 2— max(  fe, — ( fc-h  1 )) 

_gj(ci(fe+5)+5d2(fc+i)®)  ^  gi((<ii+C2(fc+|)+2a)(i7i+i)+i*(m+i)*)^ 

m=inax(fc,—  (fc+1)) 

(71) 

for —floor  (^)  <  A;  <  floor 
Special  Case  ds  =  0 


'  e>c»(*'+5)(JV  -  1  -  2max(A;, -(A:  +  1)))  if di  +  C2(fe  + 1)  +  2a  =  0 mod  tt, 
A"  rA:  a)  -  i  ^JCl(fc+i)pj[i(<il+C2(fc+i)+2a)(JV-l)]  { ”°[f(<^l+t^(fc+i)+2Q)(^-l-2inax(fc,-(fc+l)))l 

s  er  2  e-  2  ^  8in[i(di+c2(fc+i)+2Q)] 

otherwise. 

(72) 


Thus  \AlJy{k,  a)|  is  concentrated  along  the  line  di  +  C2{k  -j-  H-  2a  —  0  and 
symmetrically  spread  about  it  with  respect  to  a.  For  a  fixed  k,  the  amount  of 
spread  with  respect  to  a  can  be  measured  by 
1/{N  -  1  -  2max(fc, -(fc4- 1))). 


The  General  Case 

Although  it  seems  difficult  to  evaluate  oi)  in  the  general  case,  by 

writing 


Ai^y{k,a)  =  e^(‘^i(''+f)+i<i2(''+§)")eJ((*+2*('®^)+‘="(^+2)+2“)(V)) 

^  ^^gj((«ii+<i2(-^)+C2(fc+i)+2a))(m-(-^))+5d2(m-(-^^))*)^ 

m=max{kf—  (A:+l)) 

(73) 
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we  can  see  that  is  symmetrically  spread  about  the  line 

di  +  ^2  +  C2{k  + +  2q!  =  0  with  respect  to  ot.  For  small  d2, 

Kvik,a)]  is  concentrated  along  the  same  line.  For  a  fixed  A:,  the  amount  of 
spread  with  respect  to  a  will  depend  on  both  (iV  -  1  -  2  max(A;,  -(A:  +  1))) 
andd2. 


4.3.3  Type-Ill  Cross  Ambiguity  Function 

Recall  that  the  type-III  cross  ambiguity  function  (k,  a)  is  defined  in 
terms  of  the  type-I  and  type-II  cross  ambiguity  functions  as 


(74) 


Ai  {kl2,a) 
Aili{k-l)/2,a) 


for  even 
for  odd 


k, 

k. 


By  writing 


4f/(A;,a)  =  e^(‘=>(|)+f‘fe(|)^)e^((*+i<i2(^)+c2(|)+2«)(^)) 

m=0 

(75) 

we  can  see  that  \  A^J^^{k,  a)|  is  symmetrically  spread  about  the  line 
di  +  (h  +  ^2(2)  +  2a;  =  0  with  respect  to  a.  For  small  da, 

\-^x,y  (*^>  <^)|  is  concentrated  along  the  same  line.  For  a  fixed  k,  the  amount  of 
spread  with  respect  to  a  will  depend  on  both  {N  -  max(A:,  -k))  and  d2. 

The  Special  Case  da  =0 


Ai^,^ik,a)  =  I 
(76) 


)(JV  _  inax(A;,  -k))  if  di  +  C2(|)  +  2a  =  0  mod  tt, 

efci(|)e4[2('i»+=2(|)+2a)(W-l)]  ( sin[| (di  +C2(| )+2a)(N-max{k,-k))] \ 

y  sm[i(di-fc2(*H2Q)]  J 


Thus  |i^^^(A:,a)|  is  concentrated  along  the  line  di  +  C2(|)  -I-  2a  =  0  and 
symmetrically  spread  about  it  with  respect  to  a.  For  a  fixed  k,  the  amount  of 
spread  with  respect  to  a  can  be  measured  by  l/(Ar  -  max(A;,  -A:)). 
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5.  Detecting  Targets  in  Ci utter  using  Smoothed 

Wigner  Distributions  _ 

When  the  Wigner  Distribution  (WD)  of  a  radar  signal,  collected  by  a  pulse  Doppler 
High  Frequency  Surface  Wave  Radar  during  the  known  presence  of  a  target,  are 
displayed,  it  is  virtually  impossible  to  ‘see’  the  target  as  a  distinct  Une  in  the 
time-frequency  plane.  This  was  due  to  the  fact  that  the  strength  of  the  background 
against  which  the  target  had  to  be  detected  was  much  larger  than  the  strength  of  the 
target,  and  the  sidelobes  of  the  WD  of  the  background  and  those  of  the  interaction 
between  the  background  and  the  target  were  obscuring  the  target. 

In  general,  the  background  against  which  targets  must  be  detected  by  a  radar  will 
constitute  clutter  and  noise.  In  the  case  of  High  Frequency  Surface  Wave  Radar 
(HFSWR)  used  over  the  ocean,  the  clutter  is  due  to  Bragg  scattering  from  the  surface  of 
the  ocean  [9].  This  clutter  can  be  adequately  modeled  for  our  purposes  by  two  targets 
that  are  each  moving  with  constant  radial  velocity.  One  target  moves  towards  the  radar 
and  the  other  moves  away  from  the  radar,  with  the  magnitudes  of  radial  velocities  being 
equal.  These  are  called  the  approaching  and  receding  Braggs.  The  real  radar  signal  can 
thus  be  modeled  by  three  target  signals  -  two  of  which  represent  the  Braggs  -  plus  noise. 

Recall  that  the  WD  of  a  sum  of  two  signals  contains  cross-terms  due  to  the  interaction 
between  the  two  signals  [2].  In  the  case  of  two  chirp  signals,  the  cross-terms  are  of  an 
oscillatory  nature  and  located  midway  between  the  lines  along  which  the  auto-terms  are 
concentrated.  The  strength  of  the  cross-terms  is  proportional  to  the  square  root  of  the 
product  of  the  strengths  of  the  auto  terms. 

Approximately  speaking,  the  WD  of  a  real  radar  signal  contains  three  auto-terms  due  to 
the  Braggs  and  the  target,  three  cross-terms  due  to  the  interactions  between  the  Braggs 
themselves  and  the  Braggs  and  the  target,  and  the  interactions  of  the  Braggs  and  the 
target  with  the  noise.  In  the  vicinity  of  the  auto-term  due  to  the  target,  the  contribution 
of  the  other  terms  are  noise-like  with  amplitudes  larger  than  that  of  the  target’s 
auto-term,  thereby  obscuring  the  target. 

The  target’s  auto-term  can  be  made  clearly  visible  as  a  Une  by  smoothing  the  WD.  This 
smoothing  can  be  realized  as  a  two-dimensional  convolution  operation.  The  problem  is 
to  choose  an  appropriate  convolution  kernel  function.  In  the  following  sub-section,  we 
present  the  heuristics  of  the  derivation  of  a  smoothing  kernel  that  was  successfully  used 
on  the  real  radar  data. 

The  smoothed  WD  shows  the  target  as  a  distinct  Une  that  can  be  detected  via  a  generic 
Une  detection  method.  We  propose  taking  sUces  at  fixed  times  and  doing  a  peak 
detection  first.  Once  the  target  has  been  detected  and  a  coarse  estimation  of  its  velocity 
made  on  at  least  two  time  slices,  the  results  can  be  used  to  obtain  the  initial  velocity  and 
acceleration  parameters  of  the  target.  The  scope  of  this  study  was  Umited  to  finding  an 
appropriate  smoothing  kernel.  The  software  developed  provides  a  good  testbed  for 
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further  experimentation  with  different  choices  of  kernel  functions. 

In  the  case  of  a  Line-of-Sight  (LOS)  radar  used  over  land,  the  clutter  is  due  to 
reflections  from  the  surface  of  the  earth.  This  clutter  can  be  adequately  modeled  for  our 
purposes  by  a  stationary  target.  The  general  principles  outlined  in  this  report  for 
HFSWR  can  also  be  applied  for  an  HFLOS  (High  Frequency  Line-of-Sight)  radar. 


5.1  Heuristics  of  the  Derivation  of  the  Smoothing  Kernai 

To  choose  an  appropriate  kernel  function  we  use  our  knowledge  of  the  auto-  and 
cross-ambiguity  functions  of  chirp  signals  as  described  in  the  previous  section  and  the 
fact  that  the  ambiguity  function  is  the  two-dimensional  Fourier  transform  of  the  WD. 
Since  convolution  of  the  WD  with  a  kernel  function  is  equivalent  to  multiplication  of 
the  AF  by  the  two-dimensional  Fourier  transform  of  the  kernel  function,  we  can  choose 
an  appropriate  kernel  by  studying  the  auto-  and  cross-ambiguity  functions  of  chirp 
signals  [8]. 

We  start  by  making  some  observations  on  auto-  and  cross-ambiguity  functions  of  chirp 
signals.  Recall  that  the  AF  is  a  function  of  the  variables  time-offset  and 
frequency-offset.  In  the  following,  when  we  refer  to  the  slope  and  intersect  of  a  line  we 
associate  time-offset  with  x-axis  and  normalized  frequency-offset  with  y-axis. 

•  The  auto-ambiguity  function  of  a  chirp  signal  is  concentrated  along  a  line  through 
the  origin,  with  slope  =  -0.5  x  normalized  acceleration. 

•  The  cross-ambiguity  function  of  two  chirp  signals  with  the  same  acceleration  is 
concentrated  along  the  line  with  intercept  =  -0.5  x  difference  in  the  normalized 
initial  velocities,  slope  =  -0.5  x  normalized  acceleration. 

•  The  cross-ambiguity  function  of  two  chirp  signals  with  differing  accelerations  is 
spread  about  the  line  with  intercept  =  -0.5  x  difference  in  the  normalized  median 
velocities,  slope  =  -0.5  x  average  of  the  normalized  accelerations.  The  spread  is 
symmetric  in  the  frequency  offset  dimension  about  the  above  line.  For  small 
enough  difference  in  the  acceleration,  the  cross-ambiguity  function  is  concentrated 
along  the  above  line.  For  large  difference  in  the  acceleration,  the  cross  ambiguity 
function  gets  smeared  in  the  frequency  offset  dimension  much  like  the  Fourier 
power  spectrum  of  an  accelerating  target  gets  smeared  when  the  acceleration  is 
large. 

Based  on  the  above  we  can  make  the  following  observations  on  the  AF  of  an  HFSWR 
signal  containing  Braggs  and  one  target: 

•  For  targets  with  (magnitude)  acceleration  less  than 

Bragg  velocity/  NumPulsesPerCPI,  the  regions  where  the  cross-terms  between  the 
Braggs  are  concentrated  will  not  overlap  with  the  region  where  the  target’s 
auto-term  is  concentrated. 
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•  For  targets  with  (magnitude)  median  velocity  greater  than  3  x  Bragg  velocity,  the 
regions  where  the  cross-terms  between  the  Braggs  and  the  target  are  concentrated 
will  not  overlap  with  the  region  where  the  target’s  auto-term  is  concentrated. 

Therefore,  for  targets  satisfying  the  above  conditions  on  the  acceleration  and  the 
median  velocity,  we  can  suppress  all  cross-terms  while  preserving  the  target’s  auto-term 
simply  by  setting  the  value  of  the  AF  to  zero  over  the  region 

frequency-offset  >  Bragg  velocity  - 

0.5  X  Bragg  cross  term  region’s  width  in  the  frequency-offset  dimension. 

After  suppressing  the  cross-terms  on  the  AF,  we  can  take  the  two-dimensional  Fourier 
transform  to  obtain  a  smoothed  WD.  A  more  detailed  description  and  discussion  of 
different  methods  are  given  in  [8]. 

5.2  Analysis  of  Real  Radar  Data 

Real  HF  radar  data  was  analysed  using  the  Fourier  transform  method  and  the  type-HI 
Wigner  distribution  method.  The  data  had  been  collected  by  a  Pulse  Doppler  High 
Frequency  Surface  Wave  Radar  that  uses  a  10-element  Unear  receiving  antenna  array. 
The  data  had  been  collected  during  the  known  presence  of  a  target.  The  radar  carrier 
frequency  is  5.672  MHz  and  the  pulse  repetition  frequency  is  9.17762  Hz.  There  are 
ten  trails  and  each  trial  corresponds  to  a  block  of  256  pulses.  The  detailed  description 
of  the  radar  is  given  by  Chan  [9,10]. 

The  results  of  our  analysis  are  presented  in  Figures  1-20.  The  results  for  the  type-HI 
Wigner  distribution  of  the  time  series  containing  the  aircraft  are  compared 
simultaneously  with  the  Fourier  transform  of  the  aircraft  signal.  The  results  show  that 
whenever  the  target  was  detectable  by  the  Fourier  method,  that  is  in  trials  3, 4  and  8,  the 
target  was  detectable  also  by  the  smoothed  type-HI  Wigner  distribution  method.  For  the 
trials  3,  4  and  8,  the  aircraft  is  moving  with  a  constant  radial  velocity.  In  the  other  trials, 
that  is  in  trails  1, 2  and  5,  the  target  was  not  detectable  by  the  Fourier  method  but  the 
target  was  detectable  by  the  smoothed  type-HI  Wigner  distribution  method  to  varying 
degrees  of  success.  For  the  trials  1 ,  2  and  5,  the  aircraft  is  moving  with  a  constant  radial 
acceleration.  Trails  1,  2  and  5  show  that  there  are  limitations  and  shortcomings  to  the 
Fourier  transform  method  to  detect  accelerating  targets  because  of  the  phenomenon 
known  as  Doppler  smearing.  On  the  other  hand,  the  time-frequency  analysis  has  the 
capabiUty  to  detect  the  target,  since  its  resolution  abiUties  are  not  reduced  in  the 
presence  of  a  Doppler  signature. 

An  important  contribution  of  the  present  work  is  the  use  of  the  type-HI  WD  rather  than 
the  type-I  WD  which  has  been  used  by  many  other  researchers.  When  the  type-I  WD  is 
used,  the  range  of  unambiguously  measurable  normalized  velocities  is  tt.  Moreover, 
targets  that  are  tt  radians  away  from  the  clutter  region  in  the  spectral  domain  will  get 
masked  by  the  clutter  and  cannot  be  detected.  The  type-HI  WD  helps  us  to  overcome 
these  problems.  When  the  type-HI  WD  is  used,  the  range  of  unambiguously  measurable 
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normalized  velocities  is  27r,  and  target  that  are  tt  radians  away  from  the  clutter  region 
can  be  detected. 
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Figure  15:  The  Fourier  transform  of  the  aircraft  signal. 


Discrete  Time 


Figure  16:  The  type-ill  Wigner  dstrlbutiort  of  the  time  series  containing  the  aircraft. 
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the  aircraft. 


6.  Conclusion 


This  report  examined  the  possibiUty  of  improving  detection  of  low-altitude  aircraft 
using  High  Frequency  Surface  Wave  Radar.  Because  we  are  dealing  with  high  speeds 
and  accelerations,  the  time-frequency  analysis  has  the  capabiUty  to  improve  the 
detection,  since  its  resolution  abilities  are  not  reduced  in  the  presence  of  a  spread 
Doppler  signature. 

The  method  involved  use  of  the  Wigner  distribution,  together  with  a  ambiguity  function. 
The  results  demonstrated  superior  performance  for  all  cases  considered.  Although 
Wigner  distribution  has  for  a  long  time  been  known  to  possess  interesting  properties,  its 
success  in  solving  actual  signal  analysis  problems  has  been  limited.  The  Wigner 
distribution  has  been  generally  used  as  a  signal  interpretation  tool  only.  The  study  has 
employed  the  Wigner  distribution  for  solving  a  particularly  challenging  problem, 
demonstrating  that  use  of  time-frequency  analysis  can  improve  detection  performance. 

The  results  show  that  whenever  the  target  was  detectable  by  the  Fourier  transform 
method,  the  target  was  detectable  also  by  the  smoothed  type-IH  Wigner  distribution 
method.  In  the  other  trials  the  target  was  not  detectable  by  the  Fourier  transform 
method  but  the  target  was  detectable  by  the  smoothed  type-in  Wigner  distribution 
method  to  varying  degrees  of  success. 

Another  important  contribution  of  the  present  work  is  the  use  of  the  type-M  Wigner 
distribution  rather  than  the  type-I  Wigner  distribution  which  has  been  used  by  many 
other  researchers.  When  the  type-I  Wigner  distribution  is  used,  the  range  of 
unambiguously  measurable  normalized  velocities  is  tt.  Moreover,  targets  that  are  tt 
radians  away  from  the  clutter  region  in  the  spectral  domain  will  get  masked  by  the 
clutter  and  cannot  be  detected.  The  type-IH  Wigner  distribution  helps  us  to  overcome 
these  problems.  When  the  type-IH  Wagner  distribution  is  used,  the  range  of 
unambiguously  measurable  normalized  velocities  is  2n,  and  target  that  are  tt  radians 
away  from  the  clutter  region  can  be  detected. 
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Annex  A 

The  Wigner  Distributions  of  a  Discrete-Time  Chirp 
Signal 


In  this  appendix  we  derive  the  type-I,  type-H,  and  type- III  Wigner  distributions  for  the 
discrete-time  chirp  signal 


(A.1) 


f  ^'(6in+|62n2) 

I  0 


ifO<n<  (AT-l) 
Otherwise. 


A.1  Type-1  Wigner  Distribution 

Recall  that  the  type-I  Wigner  distribution  W/ (n,  0)  is  defined  as 


(A.2)  Wg  (n,  6)  =  ^  s{n  -1-  k)s*{n  —  k)e 

it 

where  n  is  integer- valued  and  9  is  real- valued. 

It  is  easy  to  see  that  the  signal  product  term  s(n  -|-  k)s*(n  -  k)  is  zero  outside  the 
ranges  0  <  n  <  (JV  -  1)  and  -  imn(n,  N-l-n)<k<  min(n,  iV  -  1  -  n).  Thus 
Wg(n,  0)  =  0  outside  0  <  n  <  (AT  —  1). 

Denote 


I  =  min(n,  JV— 1— n)  =  < 
(A.3) 


n 

{N-l)/2 
N  —  1—n 


ifn<iN-l)/2, 

if  N  is  odd  and  n  =  {N  —  l)/2, 

if  n  >  (AT  —  l)/2. 


In  terms  of  this,  for  0  <  n  <  (Af  -  1), 


i 

(A,4)  W^{n,9)=  e^(«>i("+*)+i»2(n+fc)")e-i(fri(n-fc)-(-§ft2(n-fe)2)g-,-2W^ 

ik=-i 

I 

(A.5)  =  ^  fj(tn2k+ib24nk)^-j2k6^ 

k=-l 

I 

(A.6)  =  ^  g-3{9-{bi+b2n))2k ^ 

k=-l 
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which  is  a  sum  of  a  geometric  series  that  can  be  easily  evaluated. 
We  first  observe  that  the  sum  can  be  written  as 


(A.7) 

k=-l 

where  a  =  2{6  —  {bi  +  i>2«))  and  then  evaluate  it  to  be 


(A.8) 


21  +  1 
sinfa(;+|)] 
sm(f) 


if  a  =  0  mod  27r, 
otherwise. 


Thus,forO<n<  (iV-1), 


Wjin,9) 

(A.9) 


2  min(n,  JV  —  1  —  n)  +  l  if0  =  6i+  b2n  mod  tt, 

sin[(g-(fci+&2n))(2mm(n,W-l-n)+l)]  otherwise. 

8in[^-(&i+&2Ti)] 


A.2  Type-1 1  Wigner  Distribution 

Recall  that  the  type-II  Wigner  distribution  Ws^(n,  9)  is  defined  as 


(A.10)  Wj\n,9)  =  ^s(n  +  A:  +  l)a*(n  - 

k 

where  n  is  integer- valued  and  9  is  real- valued. 

It  is  easy  to  see  that  the  signal  product  term  s(n  +  k  +  1)5* (n  —  k)  zero  outside  the 
ranges  0  <  n  <  (iV  -  2)  and  -  inin(n  +  l,N-l-n)<k<  inm(n,iV  -  2  -  n). 
Thus  W/^(n,  e)  =  0  outside  0  <  n  <  (iV  -  2). 

Denote 


I  =  nain(n,  JV— 2—71) 


(A.11) 


n 

<  (N-2)/2 
N-2-n 


ifn<(N-2)/2, 

if  N  is  even  and  n  =  {N  —  2)/2, 

ifn>  (iV-2)/2. 


In  terms  of  this,  for  0  <  n  <  (JV  —  2), 
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w^\n,e) 


(A.20) 


2  min(n,  JV  —  2  —  n)  +  2  ii9  =  hi  +  ^b2  +  b2n  mod  2n, 

—(2 inm(Ti,  iV  —  2  —  n)  +  2)  if  0  —  ir  +  bi  +  ^b2  +  h^n  mod  27r, 

sin[(g-(ti+f62+&2n))(2min(n,jy-2-n)+2jl  otherwise, 
sin  1^— (6i 62+?>2T^)] 


It  is  useful  to  write  this  as 


wl\n,e)  = 


(A.21) 


2min(n  +  ^,  JV  —  1  —  (n  +  ^))  +  1  if  ^  =  6i  +  h2{n  +  ^)  mod  27r, 

—  2min(n  +  5,  JV  —  1  —  (n  +  5))  +  ll  if  5  =  tt +  61 +  62(^1  + 1) 

™  (g-(fel+fr2(n+|)))(2min(n+^,JV'-l-(n+-;))+l)]  otherwise. 

mn  [5— (61+62(04-^))] 


A.3  Type-Ill  Wigner  Distribution 

Recall  that  the  type-III  Wigner  distribution  Wg^^{n,  0)  is  defined  in  terms  of  the  type-I 
and  type-n  Wgner  distributions  as 


(A.22) 


Wg  (n/2, 0)  for  even  n, 

((n  -  1) /2, 0)  for  odd  n. 


For  n even  and  0  <  n  <  2(JV  —  1), 


f  2min(f,JV-l-f)  +  l  if0  =  6i  + 

'^l^\n,0)  =  l  8in[(g-(6i+^&2n))(2min(f,JV-i-f)+i)]  otherwise. 

(A.23) 

For  n  odd  and  0  <  n  <  2(JV  —  1), 


if  0  =  61  +  hhn  mod  tt, 


Wj^\n,0)  = 


(A.24) 


2min(f,JV-l-f)  +  l 
-  ■2min(s,Ar-l-S)  +  l] 

siii[(g-(&i4-j;>2n))(2inm(f,Ar-l-f)+l) 

sm[d-(6i+i62T>)] 


if  0  =  61  +  ^b2n  mod  27r, 
if0  =  n  +  bi  +  ^b2n  mod  27r, 

otherwise. 


Combining  the  above,  for  all  0  <  n  <  2(JV  —  1), 
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(A.25) 


2min(|,iV-l- 2)^-1 
(-1)"  [2min(f,Ar-l-a)  +  i] 
»in[(^-(&i+l62n))(2min(°,W--l-a)+l)] 
8in[d-(fri+i62n)] 


if  0  =  6i  +  ^b2n  mod  27r, 
if  0  =  TT  +  6i  +  56271  mod  277, 

otherwise. 
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